Stimulated emission of two photons in parametric amplification 
and its interpretation as multi-photon interference 
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Stimulated emission of two photons is observed experimentally in the parametric amplification 
process and is compared to a three-photon interference scheme. We find that the underlying physics 
of stimulated emission is simply the constructive interference due to photon indistinguishability. So 
the observed signal enhancement upon the input of photons is a result of multi-photon interference 
of the input photons and the otherwise spontaneously emitted photon from the amplifier. 
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Stimulated emission, first proposed by Einstein [l[ to 
explain the blackbody radiation spectrum, is the main 
process in laser operation. It provides the optical gain of 
an active medium and is responsible for the coherence of 
laser light [2j. Although the process was studied exten- 
sively as an amplification process of a classical wave field 
as early as in 1955 fj|, its effect on the nonclassical state 
of light was only investigated not long ago [4j .especially 
in the contest of quantum state cloning [5|, [y, 0| ■ 

Fundamentally, stimulated emission occurs at single- 
photon level, i.e., it is seen as the creation of an identi- 
cal photon to an incoming photon. However, the same 
photon can also be produced even without the incoming 
photon, due to spontaneous emission. Thus, the exis- 
tence of the input photon will enhance the production 
rate, as compared to the case without the input photon. 
Indeed, in recent study of stimulated emission by single 
photons, a doubled rate is observed in photonproduction 
that is correlated to the input photon [1, 0, 0] . But the 
above picture is only phenomenological and does not tell 
the underlying physics. So what fundamental physical 
principle governs this phenomenon? 

If we make a detailed analysis of the enhancement due 
to stimulated emission, we find that it bears some resem- 
blance to the famous Hanbury Brown- Twiss effect [9j, 
i.e., the photon bunching effect of thermal light pj. llOj]: 
the enhancement factor in both cases is one-fold and the 
temporal profile is the same. The photon bunching ef- 
fect, as pointed out by Glauber [ll| in 1965, is in essence 
a two-photon interference effect. This suggests that the 
underlying physics in stimulated emission is simply multi- 
photon quantum interference. Recently, Khan and How- 
ell [HJ and Irvine et al. [l3[ utilized a beam splitter 
and the two-photon Hong-Ou-Mandel interference effect 
[Ti . [H[ to emulate the photon cloning process observed 
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FIG. 1: Comparison between (a) the stimulated emission and 
(b) the multi-photon interference. 



in stimulated emission [6J] . This further demonstrates the 
connection between stimulated emission and the multi- 
photon interference. 

In this letter, we wish to report on an experiment in 
support of the above claim. In the experiment, we in- 
ject a two-photon state into a parametric amplifier and 
observe an enhancement by a factor of nearly two in the 
photon production that is correlated to the input pho- 
tons. Furthermore, we mimic the same phenomenon with 
a beam splitter, in a similar manor to Rcf. 15j for the 
single-photon input case. These two phenomena can be 
viewed as a generalized three-photon bunching effect and 
are a result of three-photon constructive interference due 
to photon indistinguishability. 

To understand the connection between the stimulated 
emission and multi-photon interference, we consider the 
two situations in Fig.l. The process of stimulated emis- 
sion of an iV-photon state is shown in Fig. 1(a), where N 
photons interact with an atom in an excited state. The 
atom will emit one photon regardless of the input. Total 
photon number is N + 1. Assume that the spontaneous 
emission rate is R into the same mode of the input pho- 
tons. It is known that the emission rate stimulated by 
one photon is the same. Then since each input photon 
may stimulate the excited atom, the total rate is then 
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(N + 1)R. The extra R is from spontaneous emission. 
There is an enhancement factor of N in the photon pro- 
duction rate. The case of N = 1 was observed in Ref.0]. 

In multi- photon interference in Fig. 1(b), on the other 
hand, a single photon and N photons are combined by a 
50:50 beam splitter. The probability of detecting all N+l 
photons in one side is easily calculated to be (N+1)/2 N+1 
(see below). However, when the single photon is dis- 
tinguishable from the N photons and no multi-photon 
interference occurs, we find the detection probability is 
simply 1/2 N+1 . The N+l factor is a result of con- 
structive interference of N + 1 possibilities in detecting 
the N+l photons. Each possibility corresponds to the 
situation when the single input photon is detected by a 
specific detector. (Fig.l shows two such possibilities in 
the N + 1-photon detection.) We add the amplitudes of 
the N + l possibilities before taking the absolute value 
for the indistinguishable case but we add the absolute 
values of the amplitude of each possibility for the distin- 
guishable case. The ratio between the two cases is then 
N+l. The case of N = 1 is the Hong-Ou-Mandel photon 
bunching effect [l5[ ■ A similar interference effect was ob- 
served by Ou et al. [IB] with a |2,2) input state. Notice 
that the enhancement factor here is the same as the stim- 
ulated emission in Fig. 1(a). Therefore, the spontaneous 
emission rate R corresponds to the situation when the 
input N photons to the atom are distinguishable from 
the emitted photon by the atom, so that the atom is not 
influenced by the input photons and only emits spon- 
taneously. This case is exactly the same as the case in 
Fig. 1(b) but when the single photon is distinguishable 
from the N photons. When the input N photons are 
indistinguishable from the emitted photon by the atom, 
constructive multi-photon interference leads to a factor 
of N + 1 enhancement in photon detection rate. Thus, 
the underlying physics in the stimulated emission is the 
photon indistinguishability that results in multi-photon 
interference. 

In the following, we will use parametric amplifier 
to study the stimulated emission by multiple photons. 
Mathematically, any phase insensitive linear amplifier 
can be modelled as a single mode parametric amplifier, 
which is described quantum mechanically by [l7| 
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where cjj corresponds to the internal modes of the ampli- 
fier and is the idler mode for the parametric amplifier. It 
is usually independent of a s and is in vacuum. To pre- 
serve the commutation relation, we need \G\ 2 — \g\ 2 = 1. 

At microscopic level of atoms, we have a small value 
of \g\ « 1 or \G\ ~ 1. The unitary evolution operator 
for Eq.|T|) then has the form of 

U « 1 + (gata\ + h.c.) (2) 
With a vacuum input of |0), we have the output state 

|4>^ t = #|0)«|0>+ fl |l.>®|l,-). (3) 
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FIG. 2: Schematics for studying the stimulated emission of an 
input of iV-photon state with parametric amplification. SMF: 
single-mode fiber; IF: interference filter; H,V: horizontal and 
vertical polarizations; Th: adjustable delay. 



The last term gives the spontaneous emission with a 
probability of \g\ 2 . When the input is a single-photon 
state |l s ) <£> |0j), we have 
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The probability for the emission from the amplifier is 
then 2\g\ 2 . The extra emission probability of \g\ 2 is usu- 
ally attributed to the stimulated emission, which is sim- 
ilar to the photon bunching effect [1, [lO, EH ■ 

When the input state is a two-photon state of |2 s )|0j), 
we have the output state as 
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The photon emission rate from the amplifier is now three 
times the spontaneous rate. In fact, it is straightforward 
to find that, with an ./V-photon state as the input, the rate 
of photon emission from the amplifier is N + 1 times the 
spontaneous emission rate. As the single-photon input 
case, each fold of increase in the rate can be attributed 
to the stimulated emission from one individual photon in 
the input iV-photon state. 

Notice that when the input photons are not in the same 
mode as the amplifier and thus are distinguishable from 
the photon emitted by the amplifier, the output state 
becomes 
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where N > 1. So the photon production rate is exactly 
the same as the spontaneous emission. 

The above analysis with the parametric amplifier con- 
firms the previous results obtained from the pictorial ar- 
gument based on Fig.l. Next, we will demonstrate exper- 
imentally the enhancement effect for a two-photon input 
and compare it with a three-photon interference scheme 
with a beam splitter. 

The experimental arrangement for studying the stim- 
ulated emission is sketched in Fig. 2. A mode-locked 
Tksapphire laser operating at 780 nm is frequency dou- 
bled and the harmonic field serves as the pump field for 
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FIG. 3: (a) Four-photon coincidence of ABCD detectors in 
100s and (b) three-photon coincidence of ABD detectors in 
10s as a function of the delay cTh- 



a parametric amplifier made of a 1-mm long (3— Barium 
Borate (BBO) crystal. The crystal is so oriented that it 
satisfies the type-II phase matching condition and beam- 
like fields are generated [l8j]. A small portion is split 
from the Thsapphire laser and serves as the input field to 
the signal port of the parametric amplifier. The injected 
coherent field is heavily attenuated down to a rate much 
less than one photon per pulse. But even so, the coherent 
state consists of vacuum, one-photon state, two-photon 
state, and more. So the output state is a superposition 
of the states in Eqs.([3][n]). Therefore, in order to observe 
the enhancement effect in stimulated emission by a spe- 
cific number of photons, we need to make a projection 
measurement to the corresponding states in Eqs.(|31 [5J. 
For example, for a two-photon state input, the projec- 
tion is to the second term in Eq. ([5]) . This is achieved by 
a four-photon coincidence measurement, as depicted in 
Fig. 2. Joint measurement with the idler photon is nec- 
essary to discriminate against the three-photon contri- 
bution directly from the injected coherent field. Photon 
(in)distinguishability between the input photons and the 
photon emitted from parametric down-conversion is real- 
ized by an adjustable delay Th on the coherent injection 
field. A single-mode fiber (SMF) is used to collect the 
signal field from the amplifier, in order to ensure a good 
spatial mode match. Interference filters of bandpass of 3 
nm are used for temporal mode cleaning. 

The conditions for the situations in Eqs.([5l[6|) are sat- 
isfied by adjusting the delay Tjj- When the delay is right, 
the injected coherent field pulse arrives in time with the 
pump pulse to the amplifier and the photon emitted by 
the amplifier is indistinguishable from the incoming pho- 
tons in the coherent state. But when the delay is either 
too large or too small, there is no overlap between the 
coherent pulse and the pump pulse. This is the situa- 



tion described in Eq.©. Therefore, as we scan the delay 
Th, the four-photon coincidence of A, B, C, D detectors 
should exhibit a bunching effect with a peak-to-wing ra- 
tio close to three. Fig. 3(a) shows the result of the mea- 
surement. The error bars are the statistical errors of one 
standard deviation. The solid curve is a least square fit 
of the data to a Gaussian of the form 



F(T H ) = A 1 
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where To is the center position of the peak and T c is 
related to the width of the peak. We obtain V2 = 1.81 ± 
0.15 as the enhancement factor for the data in Fig. 3(a), 
which gives 2.81 as the ratio between the peak and the 
wing. This value is close to the ideal value of three in 
Eq.© for the stimulated emission by two photons. 

In the meantime, three-photon coincidences of ABD 
detectors are also registered and are shown in Fig. 3(b). 
This measurement corresponds to the second term in 
Eq.(@} and gives the stimulated emission by one input 
photon. The Gaussian fit gives an enhancement factor of 
Vi = 0.88 ± 0.14. The peak-to-wing ratio of 1.88 is close 
to the ideal value of two in Eq.©. 

The reason for the imperfection in the experiment 
is due to mode mismatch between the input field and 
the amplifier mode, i.e., mismatch between |JV) and the 
mode for which the operator a s represents. Although 
the spatial mode is matched by the single-mode fiber 
(SMF in Fig. 2), the temporal mode is hard to match 
because the temporal coherence of the parametric down- 
conversion process is very complicated and the fields are 
not transform-limited even if the pump field is. Neverthe- 
less, we use interference filters to clean up the temporal 
profile. The full widths of the peaks in Fig. 3 is approx- 
imately 2T C = 660ps, close to the coherence time of the 
interference filters (IF in Fig. 2) of width 3 nm. 

Next we consider the situation depicted in Fig. 1(b) 
where an iV-photon state is superposed with a single- 
photon state by a 50:50 beam splitter. The output state 
for the beam splitter is given by [lj| 
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where we only write down the state for which all the 
N + 1 photons exit at one port (port 1) of the beam 
splitter. On the other hand, if the input N photons are 
distinguishable from the single photon from the other 
input port, they behave like classical particles and follow 
the Bernoulli distribution. The output state becomes 
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Therefore, the rate of detecting N + l photons in port 1 is 
N times bigger when the photons arc all indistinguishable 
than when the N photons are distinguishable from the 
one photon. As discussed before, this increase stems from 
a constructive N + 1-photon interference. 
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FIG. 4: A simple scheme for interference of N photons and 
one photon with a beam splitter. Same notations as Fig. 2. 
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weak coherent field is directed to a beam splitter to com- 
bine with the signal field (H) from the parametric down- 
conversion. In order to mimic the stimulated emission 
process shown in Fig. 2, all the experimental parameters 
such as pump power, the strength of coherent field, etc. 
are the same as those in Fig. 2. We adjust the delay Tr on 
the coherent field to ensure the temporal overlap between 
the coherent field and the down-converted photon. When 
gated on the detection of the V-photon by detector D, the 
H-field of the down-conversion is in a single-photon state. 
But because of the complicated dispersion in the para- 
metric down-conversion process, the single-photon state 
is not transform-limited. Again, interference filters are 
used to clean up some of the temporal modes but impos- 
sible to achieve the perfect match. 

We record both the four-photon coincidence of ABCD 
detectors and the three-photon coincidence of ABD de- 
tectors. The former corresponds to the N = 2 case in 
Eq.©, whereas the latter to the N = 1 case. The results 
are shown in Fig. 5. The fitted curves are very similar to 
Fig.3 but with v 2 = 1.78 ± 0.14 and v x = 0.86 ± 0.09. 

As can be seen, Figs. 3 and 5 show the same result 
within the statistical errors. This confirms our claim that 
the underlying physics in stimulated emission is nothing 
but multi-photon interference. The interference effect is 
a result of indistinguishability between the input photons 
and the photon emitted by the amplifier. 
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FIG. 5: Same as Fig. 2 but with data obtained from Fig. 4. 



The experimental arrangement with a beam splitter is 
similar to Fig. 2 and is shown in Fig. 4, where the split 
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